Inter-subband resistance oscillations in crossed electric and magnetic fields 
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Quantum oscillations of nonlinear resistance are investigated in response to electric current and 
magnetic field applied perpendicular to single GaAs quantum wells with two populated subbands. 
At small magnetic fields current-induced oscillations appear as Landau-Zener transitions between 
Landau levels inside the lowest subband. Period of these oscillations is proportional to the magnetic 
field. At high magnetic fields different kind of quantum oscillations emerges with a period,which 
is independent of the magnetic field. At a fixed current the oscillations are periodic in inverse 
magnetic field with a period that is independent of the dc bias. The proposed model considers these 
oscillations as a result of spatial variations of the energy separation between two subbands induced 
by the electric current. 
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INTRODUCTION 

The magnetotransport phenomena in high-mobility 
modulation-doped semiconductor structures are com- 
monly studied with only one populated subband (Ei), 
because the electron mobility decreases with filling sec- 
ond subband (E 2 ) due to inter-subband scattering [![. 
The latter also gives rise to magneto-inter-subband oscil- 
lations (MISO) of the dissipative resistance Q. In elec- 
tron systems with two populated subbands MISO have 
maxima in magnetic fields B satisfying the relation 0- 
5]: A12 = I ■ HlOc, where A i2 = E 2 — E\ is the energy 
separation of the bottoms of the subbands, lo c — eB/m* 
is cyclotron frequency, m* is effective electron mass and 
I is a positive integer. In contrast to Shubnikov dc Haas 
(SdH) oscillations the MIS-oscillations exist at high tem- 
perature kT > huj c . An interference of these oscillations 
with phonon-induced oscillations has been reported @. 

At small quantizing magnetic fields a finite electric 
current induces several additional nonlinear phenomena. 
At low temperatures small currents decrease consider- 
ably the resistance. The dominant mechanism, inducing 
the resistance drop, is a peculiar Joule heating (quan- 
tal heating), which produces a non- uniform spectral dif- 
fusion of electrons over the quantized spectrum. The 
spectral diffusion is stabilized by inelastic processes (" in- 
elastic" mechanis m) Ml - The heating is observed and 
studied recently [8l4ll| . At higher currents electron tran- 
sitions between Landau levels occur due to an elastic 
electron scattering on impurities in the presence of elec- 
tric field 12|, |l3j. The transitions increase the resistance, 
which was observed in electron systems both with a sin- 
gle occupied subband 
occupation 17 - 2j 



stronger magnetic fields reveal another kind of current- 
induced resistance oscillations in electron systems with 



14-lfJ and with multi-subband 



In the last case an interference of 
the magneto-inter-subband quantum oscillations (MISO) 
with the current induced inter-level scattering was re- 
ported. 

Recent investigations of the nonlinear transport in 



a single band occupation[2l|. These oscillations occur 
in electric fields that are significantly smaller than the 
one required for the current-induced Landau-Zener tran- 
sitions between Landau levels 12[. The period of the 



current-induced oscillations is found to be independent 
of the magnetic field. The oscillations are considered to 
be a result of spatial variations of the electron filling fac- 
tor (electron density Sn) with the applied electric field. 

In this paper we report an observation of current- 
induced resistance oscillations of the dissipative resis- 
tance in electron systems with two populated subbands. 
Two kinds of oscillations are detected. At small magnetic 
fields we observed resistance oscillations with a period 
proportional to the magnetic field. We found that these 
oscillations are related to the current induced Landau- 
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Zcncr transitions between Landau levels 
higher magnetic fields another type of the resistance oscil- 
lations emerges with a period that is independent of the 
magnetic field. In the paper these oscillations are studied 
at high temperatures at which only MIS-oscillations are 
present. 

Despite a similarity between the current induced oscil- 
lations with the B-independent period, which are found 
in single subband systems [2l| and the oscillations re- 



ported in this paper, there is at least one distinct fea- 
ture to distinguish the two. Namely the oscillations 
in the two-subband systems occur at high temperatures 
kT 3> huj c and, therefore, the total number of the elec- 
tron states carrying the electric current (inside the energy 
interval kT) does not oscillate with the Fermi energy (in 
other words with the total electron density n). In this 
regime the SdH oscillations are damped and in single sub- 
band systems the curren-induced oscillations are absent 
211. Thus even if both kinds of observed oscillations have 



a common origin, the oscillations reported in this paper 
are not directly (simply) related to the spatial variations 
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of the electron density 5n induced by the electric current. 
Another interesting feature is the phase of these oscil- 
lations . The oscillations appears to be quasi-periodic 
with respect to the applied current but with an apparent 
7r-phase shift with respect to the zero bias. Below we 
present our findings and provide an interpretation of the 
obtained results. 

EXPERIMENTAL SETUP 

Our samples are high-mobility GaAs quantum wells 
grown by molecular beam epitaxy on semi-insulating 
(001) GaAs substrates. The width of the GaAs quan- 
tum well is 13 nm. Two AlAs/GaAs type-II superlattices 
grown on both sides of the well served as barriers, pro- 
viding a high mobility of 2D electrons inside the well at a 
high electron density [22J. Two samples were studied with 
electron density n\,2 — 8-09 xlO 15 m -2 and mobility fi± = 
121 m 2 /Vs and fi 2 '= 73 m 2 /Vs 

The studied 2D electron systems are etched in the 
shape of a Hall bar. The width and the length of the mea- 
sured part of the samples are d =50/zm and L =450/im. 
To measure the resistance we use the four probes method. 
Direct electric current Idc {dc bias) is applied simulta- 
neously with 12 Hz ac excitation I ac through the same 
current contacts (x-direction) . The longitudinal ac (dc) 
voltage V££ (V^r) is measured between potential contacts 
displaced 450^m along each side of the sample. The Hall 
voltage Vh is measured between potential contacts dis- 
placed 50fim across the electric current in y-direction. 

The current contacts are separated from the measured 
area by a distance of 500,um, which is much greater 
than the inelastic relaxation length of the 2D electrons 
Lin — (Drin) 1 / 2 ~ 1 — 5/xm. The longitudinal and Hall 
voltages were measured simultaneously, using two lockin 
amplifiers with 10 Mf2 input impedances. The potential 
contacts provided insignificant contribution to the over- 
all response due to small values of the contact resistance 
(about lkf2) and negligibly small electric current flowing 
through the contacts. 

Measurements were taken at different temperatures 
and magnetic fields in a He-3 insert inside superconduct- 
ing solenoid. Samples and a calibrated thermometer were 
mounted on a cold copper finger in vacuum. Magnetic 
field was applied perpendicular to the 2D electron layers. 

RESULTS 

Figure [T] presents the dependence of the dissipative re- 
sistance on the magnetic field at temperature T — 4.35 K. 
At this temperature kT > Tilo c and Shubnikov dc Haas os- 
cillations are suppressed at B <0.5 T. The maximums of 
the observed magneto-intersubband oscillations (MISO) 
are due to the enhancement of electron elastic scattering, 
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FIG. 1: Dependence of the resistance R xx on magnetic field 
with no dc bias applied. Sample Nl. 
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FIG. 2: a) Dependence of differential resistance R xx on mag- 
netic field and averaged density of electric current J; (b) De- 
pendence of the resistance on the current density J at fixed 
magnetic field as labeled. Index j = ±1,±2... numerates 
Landau-Zener transitions inside lowest subband, which obey 
Eq.©. T=5.1 K. Sample Nl. 

which occurs when the Landau levels in two subbands 
are lined up with each other (state P in FigfT]). At this 
condition elastic electron transitions occur between the 
subbands, increasing the total electron scattering rate 
and, thus, the resistance. Minima of the oscillations oc- 
cur when the Landau levels in one subband are between 
the levels of another subband. In this condition the elas- 
tic electron scattering between subbands is suppressed 
(state M in FigfT]) Q. 

Figure [^fa) presents differential resistance R xx at dif- 
ferent averaged density of the electric current J = 
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T=2.1 K 
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FIG. 3: Dependence of resistance R xx on magnetic field and 
current density J. Labels +A. +B and +C indicate different 
maximums induced by dc bias. T=2.I K. Sample Nf . 



Idc/ (d = 50/z) and small magnetic fields. The differential 
resistance oscillates with the dc bias. An example of the 
oscillations is shown in figure HJb) at fixed magnetic field 
B=0.12 Tesla. The dependence is a horizontal cut of the 
2D plot and is shown by the dashed line in FiglJJa). Po- 
sition of a resistance maximum j is proportional to the 
magnetic field and satisfies the following relation: 



2eE j R ( P = j ■ huj c , 



(1) 



where Ej is the electric field (mostly the Hall electric field 
in the sample) corresponding to the maximum j , is 
the cyclotron radius of electrons in the first subband (the 
lowest subband) and j = 0, 1, 2... is an integer. Eq.([T|) de- 
scribes Landau-Zener transitions between Landau levels 
in the first subband jl~2l]. 

At a higher resolution the data shows oscillations of 
the magnitude of the maximums j = ±1 with the mag- 
netic field at B > 0.1 (T). The oscillations are periodic 
in inverse magnetic field and are in-phase with the in- 
tersubband oscillations at zero dc bias (j = 0). Similar 
oscillations are observed for the minimum between j = 
and j = ±1 maximums. These oscillations are shifted by 
phase 7r with respect to the oscillations of the maximums 
j = 0, ±1. The observed oscillations appear as an inter- 
play between the dc bias induced Landau-Zener transi- 
tions between Landau levels inside the lowest subband 
and the intersubband transitions, which are periodic in 
inverse magnetic field 1/B. At higher dc biases > 1) 
the amplitude modulation with the 1/B periodicity dis- 
appears. In particular no amplitude modulation is found 
for j = ±2,3 maximums. 




FIG. 4: (a) Dependence of resistance R xx on magnetic field 
and current density J. indicating strong correlation of fea- 
tures ±A and ±C with MISO minimums and features ±B 
with MISO maximums, (b) Dependence of R xx on current 
density J at magnetic field B=0.418 T corresponding to MISO 
maximum and at magnetic field B=0.408 T corresponding to 
MISO minimum. T= 4.7K. Sample Nl. 



Figure [3] presents a typical nonlinear response at a high 
magnetic field. The response is symmetric with respect to 
applied dc bias and is shown for the positive bias. There 
are several distinct features, which appear with the dc 
bias. The features are labeled in the figure. Firstly we 
discuss evolution of the resistance with the dc bias at 
minimum of a MIS oscillation (state M in fig[T]). When 
the dc bias is applied the resistance falls down and, then, 
develops a shoulder labeled by symbol +A. The initial 
drop of the resistance is mostly due to the quantal heat- 
ing. Further increase of the dc current leads to formation 
of a maximum labeled by symbol +C. 

When the dc bias is applied to state P (see figH]), corre- 
sponding to maximum of a MIS oscillation, the resistance 
drops much more abruptly and significantly in compar- 
ison with the previous case. At low temperatures the 
resistance drop reaches zero and forms zero resistance 
state (ZDRS) [181 |23I425| . Further increase of the dc bias 
leads to the formation of a maximum labeled by symbol 
+B. 

An evolution of the discussed features with the mag- 
netic field is shown in figure HJa) . The figure demon- 
strates that the positions of all features (±A, ±B, ±C) 
are essentially independent of the magnetic field. Figure 
IHb) presents horizontal cuts of the 2D plot through a 
maximum (B=0.418 T) and a minimum (B=0.408 T) of 
the inter-subband quantum oscillations. 

Figure [5] presents an overall behavior of the quantum 
oscillations in a broad range of magnetic fields and dc 
biases. The data was obtained on sample N2. The fig- 
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FIG. 5: Evolution of differential resistance with magnetic field 
and current density in broad range of magnetic fields. White 
straight lines indicate Landau-Zener transitions which obey 
Eq.|T]). Upper panel presents horizontal cut through MISO 
maximum at B=0.548T (gray line) and cut through MISO 
minimum at B=0.532 T (black line). Sign +(-) indicates re- 
gions of current density J, inside which the current induced 
oscillations have 0(180) degree phase shift with respect to 
MIS-oscillations at J=0A/m. Right panel presents two verti- 
cal cuts of the 2D plot taken at current densities as labeled. 
Magnetic filed dependence at J=3.03 A/m indicates strong 
reduction of the resistance oscillations at B < B c inside the 
region corresponding to Landau-Zener transitions. T=5K. 
Sample N2. 

ure shows the crossover of the intraband Landau-Zener 
transitions, obeying Eq.([T]), and the oscillations marked 
as ±A, ±B, ±C, which have the MISO periodicity. The 
apparent crossover occurs near the Landau-Zener transi- 
tion corresponding to j — ±1. Namely the oscillations 
with 1/B MISO periodicity occurs at magnetic fieds B c 
corresponding to 

hco c > 2eE 1 R i c L) . (2) 

At smaller magnetic fields B < B c the oscillations are 
significantly reduced. Two vertical cuts of the 2D plot 
taken at different currents are shown in the right panel of 
Fig. The curve taken at J=3.03 A/m shows the strong 
reduction of the oscillations at B < B c in a compari- 
son with the MISO at J=0A. Thus the main intraband 
Landau-Zener transition {j = ±1) forms a boundary be- 
low which the current induced oscillations with 1/B in- 
tersubband periodicity are strongly damped. 

The upper panel of Fig. [5] shows two horizontal cuts 
of the 2D plot. Black solid line presents dependence of 
the resistance R xx on dc bias taken at B=0.532(T) corre- 
sponding to a mimumum of MISO. The grey line presents 
the dependence taken at B = 0.548 (T) corresponding to 
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FIG. 6: Positions of resistance maximums and different mag- 
netic fields and current density. Two kind of oscillations are 
observed: in magnetic fields at and below B c , which satisfy 
Eq.©) the maximums correspond to Landau-Zener transi- 
tions in lowest subband that obey Eq.fTJ. Solid straight lines 
at j — ±1, 2, and 3 represent the equation. At B > B c the re- 
sistance maximums follow the vertical solid lines representing 
features ±A, ±B and ±C shown on Fig |3l4l5l The crossover 
between two kind of oscillations occurs at B — B c presented 
by line j — ±1. Sample Nl. 

a MISO maximum. The two curves intersect at 8 points. 
These intersections marks the regions at which the oscil- 
lations with MISO periodicity changes their phase by ir. 
At the intersections the oscillations are nearly vanished. 
Sign " +" indicates the region between two intersections 
in which the oscillations are in-phase with the MISO, 
whereas sign "-" indicates the regions in which the oscil- 
lations are shifted by phase ir with respect to the MISO. 

Figure [6] presents an accurate position of the resis- 
tance maximums with 1/B periodicity at different cur- 
rents and magnetic fields for sample Nl. The figure indi- 
cates clearly that at B = B c (j = ±1) the resistance max- 
imums follow the main Landau-Zener transition j = ±1 
whereas at B > B c the maximums are nearly indepen- 
dent on magnetic field (features ±A, ±B, ±C). The solid 
lines j = ±1 mark the boundary between the two kinds 
of oscillations. The lines obey eq.([T]) at j = ±1 with the 
cyclotron radius R\ corresponding to the lowest subband. 
The complete theory of the current induced oscillations of 
the resistance of 2D electron system with two populated 
subbands is not available in a general case. The case 
of a bilayer electron system with two closely spaced and 
almost equally pop ulated electronic subbands has been 
studied recentlvfl7l [l9j . These results are in qualitative 
agreement with the present data at small magnetic fields 
B < B c . 

At high magnetic fields B > B c figures[S]and[n]presents 
new kind of the current induced quantum oscillations. A 
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FIG. 7: Dependence of resistance on inverse magnetic field at 
different dc biases as labeled. T=5 K. Sample N2. 



striking feature of these oscillations is the independence 
of the position of these oscillations on magnetic field. An 
interesting property of these oscillations is the region in 
which the oscillations occur. Figures 0311] show that these 
oscillations start at the line corresponding to Landau- 
Zener transitions at j = ±1 in the lowest subband and 
propagate to higher magnetic fields. Another interesting 
property is an apparent quasi-periodicity of the oscilla- 
tions with applied current. Namely the features ±A, ±B, 
±C are displaced by about the same value of the electric 
current density from each other: SJ ~1.27 (A/m). The 
phase of the oscillations is shifted by ir with respect to 
zero dc bias. It seems strange that the MIS-oscillations 
(J=0 (A/m) are not a part of this periodic set. 

Figure [7J demonstrates the 1/B periodicity and the 
phase of the current induced oscillations at different dc 
biases as labeled. The figure indicates that oscillations 
at J=1.97 (A/m) ( B+ feature) are in phase with MISO, 
whereas oscillations at J=0.575 A/m (A+ feature) are 
shifted by ir with respect to MISO. Figure [7] shows also 
the strong reduction of the oscillations at J=0.971 A/m. 
At this current the oscillations change phase by tt. The 
current corresponds to the intersection of two curves 
shown in upper panel of FigfSJ 

The 1/B periodicity of the oscillations and the mag- 
netic field independence of the electric current Idc, in- 
ducing the oscillations at B > B c , indicates a similarity 
of these quantum oscillations with the current induced 
quantum oscillations reported recently in Ref . [2l| . Below 
we consider a model, which is, in many respects, analo- 
gous to one described in Ref. 2l| . The model reproduces 



MODEL AND DISCUSSION 

Current induced quantum oscillations with 1/B period- 
icity were recently observed in 2D electron systems with 
a single occupied subband [HI]. The oscillations occur in 
a strong magnetic field at which Shubnikov de Haas oscil- 
lations (SdH)[26| are well developed. With respect to the 
electric current the oscillations are periodic with a period 
that is independent on the magnetic field. The proposed 
model considers the oscillations as result of a variation of 
the electron filling factor with the dc bias. In contrast to 
SdH oscillations, the variation appears across the sample 
and is related to a spatial change of the electron density 
Sn. If the change Sn is comparable with the number of 
electron states in a Landau level hq = m/iirh 2 ) ■ Hlu c , 
then one should expect a variation of the electron resis- 
tivity. The spatial variation of the resistivity leads to 



oscillations of the sample resistance 21 



MIS-oscillations occur due to a periodic enhancement 
of the inter-subband scattering, when Landau levels in 
two subbands are line up as shown in Fig|T] MISO have 
maxima in magnetic fields B satisfying the relation : 
A12 = ihuj c , where A12 = E2 — E\ is the energy separa- 
tion of the bottoms of the subbands and i is an integer. 
In contrast to SdH oscillations the MIS-oscillations ex- 
ist at high temperature kT > Huj c and are insensitive to 
variations of the Fermi energy and/or electron density n 
for non-interacting 2D carriers. 

For interacting electron systems the situation is dif- 
ferent. Recent direct experiment indicates that gap Eq 
between conducting and valence bands of 2D electron 
systems formed in GaAs quantum wells depends consid- 
erably on the electron density n 27J. This observation 



the main properties of the observed quantum oscillations. 



opens a way to consider the dependence of the energy 
separation between two subbands A12 on the electron 
density as a mechanism leading to the current induced 
quantum oscillations in magnetic fields B > B c . Indeed 
the experiment Ref. 27] demonstrated about one percent 
change of the gap Eq at a Hall voltage Vh—75 (mV) in 
magnetic field B=0.3 (T). The Hall voltage is compara- 
ble with the one observed in our experiment: Vh ~ 50 
mV at B=0.35 (T) and J=4 (A/m). At B=0.35 (T) the 
phase of the MISO 2-7rAi2/?ia; c « 2-7T • 30 requires about 
3 percent change of the inter-subband energy separation 
A12 to make an additional MIS-oscillation cycle. The 
comparison indicates a feasibility of the proposed mech- 
anism, taking into account that in our samples the GaAs 
quantum well is sandwiched between conducting layers, 
which enhance significantly the electron screening and, 
therefore, the variations of the electron density Sn with 
the dc bias [HI. 

In the model described below we assume that the 
dc bias-induced variation of the electron density Sn(r) 
changes the energy separation A 12 (n) between two sub- 
bands across samples. Since relative variations of the 
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FIG. 8: Dependence of the electric potential on position y in 
the direction perpendicular to the electric current in strong 
magnetic field. Line V^ D describes the potential in GaAs 
quantum well, in which strong Hall effect is developed. Line 
Vff L describes the potential in the highly resistive superlattice 
layer, in which the Hall voltage is negligibly small due to the 
negligibly small current in the layer. 



electron density is small Sn/n <C 1, we will consider only 
the linear term of the dependence Ai2(n): 



A 12 (n) = A + j8n(r), 



(3) 



where Ao is the energy separation at zero dc bias and the 
parameter 7 is a constant. The following consideration is 
qualitatively similar to the model described in detail in 
Ref . pH . Below we describe the main parts of the model 
omitting some details. 

The conducting 2D electron system in the GaAs quan- 
tum well is sandwiched between two layers of AlAs/GaAs 



superlattices (SL) of the second kind [22J. The param- 
eters of the superlattices are adjusted to set the system 
close to a metal-insulator transition. At this condition 
the barely-conducting SL layers efficiently screen electric 
charges but do not contribute considerably to the overall 
conductivity of the structure. Electric contacts connect 
the GaAs and the SL layers. Thus the system is consid- 
ered as a set of parallel conductors. At zero magnetic 
field the distribution of the electric potential driving the 
current is the same in all layers due to the same shape of 
the conductors. That is to say at B=0 the potential dif- 
ference between different layers is absent. In the poorly 
conducting SL layers the electric current is several order 
of magnitude smaller than the one in the highly conduct- 
ing GaAs quantum well. 

The layers have a different distribution of the electric 
potential in a strong magnetic field, at which uj c t?® J> 1 
and oj c rf r L <C 1, where and rf r L are transport times 
in the GaAs and in the SL layers. At uit^ 1 the elec- 



tric field in the GaAs layer is almost perpendicular to 
the current due to the strong Hall effect. In contrast the 
very small electric current in the SL layer induces a Hall 
voltage, which is negligible. The Hall voltages are shown 
in Fig|S](a) for small currents ( linear response). Figure 
|S{b) presents distribution of electric charges in the struc- 
ture. Electric charges are accumulated near the edges of 
the 2D highly conducting GaAs layer, inducing the Hall 
electric field Eh- The charges are partially screened by 
charges accumulated in the conducting SL layers. 

Due to the small Hall voltage E H L and the absence 
of the electric current across the system the change of 
the electric potential <p SL (y) in the SL layer is negligibly 
small. Below we consider the potential 



lSL 



as a con- 



stant. Due to a finite screening length A s in the SL layer 
the charge accumulation occurs at a distance d ~ X s . Be- 
low we approximate the charge distribution by a charged 
capacitor with an effective distance d e ff between con- 
ducting plates. 

The proposed model considers a long 2D Hall bar with 
a width L y 2&, 2t| . Electric current is in ^-direction and 



the Hall electric field is in y-direction. In a long conduc- 
tor the electric field E = (E x ,E y ) is independent on x, 
due to the uniformity of the system in x direction: 



dE x dE v 



= 



dx dx 

For a steady current Maxwell equations yield: 
8E X dE v 



(4) 



dy 



dx 



(5) 



Eq.Q and eq.© indicate, that the x component of the 
electric field is the same at any location: E x = E =const. 

Boundary conditions and the continuity equation re- 
quire that the density of the electric current in y direction 
is zero: J y = and therefore, 



E r 



PxxJx E y PyxJx 



(6) 



where p xx and p yx are longitudinal and Hall components 
of the resistivity tensor 30] . We approximate the MIS- 
oscillations of the resistivity by a simple expression [5| : 



. . ., . ,2ttAi2,, , s 

Pxx{n(y)) = p D [l + a ■ cos{— ) (7) 

nu) c 

where po is Drude resistivity, and a describes the ampli- 
tude of the quantum oscillations. 

An electrostatic evaluation of the voltage between con- 
ducting layers, shown in FigJSJb), yields: 



eSn(y)d ef f 
2ee 



(8) 
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where (j) 2D and <p SL are electric potentials of the GaAs 
(2DEG) and superlattice (SL) layers, and e is permittiv- 
ity of the SL layer. Expressing the electron density 5n 
in terms of electric potential (p 2D from Eq.© and sub- 
stituting the relation into eq. © and then into eq. © one 
can find dependence of the resistitivity on the electric 
potential: Pxx{4> 2D )- 

The relation E y = —d(j) 2D /dy together with eq.© 
yields: 



i2D 



Pxx(0 2D ) = PyxE 



dy 



(9) 



Separation of the variables (f> 2D and y and subsequent 
integration of eqU] between two sides of the 2D conduc- 
tor (y-direction) with corresponding electric potentials 
fa and fa yield the following result: 



p D (fa - fa + f[sin[&(fa - fa)} 
xcos[§(0 2 + fa) + 9 }]) = p xlJ EL y 

(3 = 4Tre ej/{ed eff n ), 
6 = 2nA /nu c - ^ SL , 



(10) 



where L y is a width of the sample. Taking into account 
that longitudinal voltage is V xx = EL X , where L x is a dis- 
tance between the potential contacts, and the Hall volt- 
age Vh = fa - fa = - / Eydy = -p yx I (see eqlSj), the 
following relation is obtained: 



Vx: 



Rd(I 



fiPx 



(11) 

, where Rd = L x pd/L v is Drude resistance. 

Eq llll is simplified further for two cases corresponding 
to a minimum and a maximum of MIS-oscillations. In 
these cases the voltage 4> 2D (5y) — <j) SL is expected to be an 
asymmetric function of the relative position 5y = y — yo 
with respect to the center of the sample yo (as shown in 
fig. [SJ and, thus, fa - 4> SL = — (fa - 4> SL ) and the argu- 
ment of the cosine in eq Jill becomes to be independent 
on the electric current. In these cases the differential 
resistance r xx = dV xx /dI is found to be 



I 27rAg 

= Rd\\ + a • cos(2tt— ) ■ cos{— )J 

1q ncu c 



(12) 



where Iq = e 3 hd e f fn/eeomj. 

Eq. rrj] demonstrates oscillations of the differential re- 
sistance with the electric current. The period of the os- 
cillations Iq does not depend on the magnetic field in 
accordance with the experiment. The amplitude of the 
MIS-oscillations is strongly modulated by the dc bias. In 
particular at I — Jo/4 the amplitude is zero. At this 
node the 1/B periodic oscillations change phase by n. 
The strong amplitude modulation with the dc bias and 
the 7r phase shift at a node agree with the experiment. 



Following from Eq. (fT2|) positions of the nodes and anti- 
nodes of the oscillations with respect to the current Idc 
do not agree with the experiment. In accordance with 
Eq. (TT2")) the nodes occurs at 



2i-l:i 



1,2,3..., ^ 

where k is a node index. Upper panel of Figj5] shows 
nodes at 0.22, 0.93, 2.41 and 3.91 A/m. Thus the relative 
positions of the nodes observed in the experiment do not 
follow the node positions (or index k) in Eq.(|13p. Below 
we show that the disagreement is reduced significantly 
taking into account the Joule's heating. 

The discussed above model does not take into account 
the dc heating of the 2D electrons. The Joule heating in 
systems with a discrete spectrum (quantal heating) has 
a peculiar form providing strong impact on the electron 
transport [§| . In electron systems with two subbands oc- 
cupied the quantal heating inverts the MIS-oscillations 
Tol 11 1 . A quantitative account of the heating will not 
be done in this paper. Instead we will use an analyti- 
cal approximation of the heating which is valid for two 
subbands with equal electron population. As shown be- 
low the approach yields the positions of the nodes which 
agree with the experiment. 

The expression for the resistivity of 2D electron sys- 
tems with two equally populated subbands in crossed 
electric and quantizing magnetic field reads [Tl| 

Pxx^PD[l + exp(-^-)^(l 



cos 



'2ttA 



~)] 



O = 2 f J l 



(14) 



where r q is quantum scattering time, Tj n and r tr are in- 
elastic and transport scattering times. To account the 
heating we replace Eq.([7]) by Eq. (fT4"]) and evaluate dif- 
ferential Eq.© numerically with fitting parameters ap- 
proximating the experimental data. Due to a quite rough 
approximation of the heating, the fitting parameters may 
deviate significantly from actual physical values. To find 
the fitting parameter corresponding to the inelastic scat- 
tering time we use that the second term of Eq. (fT4"f is zero 
at Q=l/3 11]. Assuming that at a small dc bias and 
low temperatures the quantal heating dominates 0, fll| , 
we related the first node shown in FigJS]at J=0.22 A/m 
to the condition Q=l/3. This yields r, n = 1.8 ns at 
B=0.53T. Using this value we solved Eq.© numerically. 
The result is shown in Figj9](a). At small dc bias J «0.17 
A/m the figure demonstrates the oscillation node, in- 
duced by the heating with a small contribution from the 
variation of the band separation A12. Other nodes oc- 
cur at considerably higher dc biases and are shifted with 
respect to the nodes shown in Fig[9]Jb)), which obtained 
by the numerical evaluation, ignoring the quantal heating 
(Q=0). 

At Q > 1/3 the heating not only shift the nodes but 
also inverts of the oscillations induced by variations of 
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FIG. 9: (a) Numerical simulation of the dependence of dif- 
ferential resistance on dc bias at B=0.53 T. Fitting param- 
eters used in the numerical simulation: Ti n — 1.8 ns, r 9 =2.5 
ps and rj r =45 ps; electron density n=8.09-10 15 1/m 2 ; effec- 
tive screening length d e //=30 nm; parameter 7=1-10 -37 Jm 2 
(see Eq.([3})- (b) Numerical simulation of the dependence of 
differential resistance on dc bias with the same fitting pa- 
rameters as in (a) but without dc heating: Ti„—0 ns (Q=0). 
Filled (open) circles present evolution of a MISO maximum 
(minimum) with the dc bias 



the band separation. Namely shown in FigJ^a) the max- 
imum at J=1.75 A/m is a result of the dc bias induced 
evolution of the MISO maximum at J=0A/m. Without 
the heating the MISO maximum evolves into a minimum 
at J=1.65 A/m shown in Figj9^b). Thus the heating 
inverts minimums to maximums and visa versa. The in- 
version is directly related to the sign change of the second 
term in Eq.([Hj) at Q=l/3. 

The heating and the variation of the band separation 
affect differently the maximums and minimums of MIS- 
oscillations. Conversely, quantal heating decreases the 
resistance at any magnetic field. A variation of the re- 
sistance, induced by the change of the band separation, 
depends on the magnetic field. At a maximum (state P 
in Fig.([T])), a variation of A i2 destroys the level align- 
ment decreasing the inter-band scattering and, thus, the 
resistance. At a minimum (state M in Fig.([T]), a varia- 
tion of A12 improves the level alignment and increases 
the inter-band scattering and the resistance. Thus at 
a MISO maximum both the heating and the variations 
of the band separation decreases the resistance whereas 
at a MISO minimum two mechanisms work against each 
other. In result the drop of the resistance at a MISO 
maximum is considerably stronger the one at a MISO 
minimum. In fact the shoulder (feature +A in FigOJ) is 
a result of the competition between two mechanisms at 
a MISO minimum whereas ZDRS states, developed from 



FIG. 10: Position of nodes of dc bias induced oscillations 
shown in Fig[9]and Fig[5]at different node index k. Filled tri- 
angles present nodes, which are obtained numerically without 
heating and obey Eq. (|13p (solid line). Account of the heat- 
ing (filled squares) improves significantly agreement with the 
experiment (open circles). Dashed line is a shift of the solid 
line to the right by two units (see text for detail). 



MISO maximums, is a strong indication of the joint de- 
crease of the resistance due to both mechanisms. The 
behavior is reproduced in the proposed model. Indeed 
FigEJa) shows that the initial drop of the MISO max- 
imum is considerably stronger than the decrease of the 
MISO minimum with the dc bias. 

Figure [10] presents a comparison of the positions of 
oscillation nodes, obtained in the model, with the exper- 
iment. For the purpose of a comparison, the node posi- 
tions are plotted versus the index k, which is defined in 
Eq. flTB")) . Without the heating nodes of oscillations obey 
Eq|T3l Filled triangles demonstrate this behavior. When 
the heating is on ( filled squares) the first node ( k=l) is 
due mostly to the heating. The following nodes ( k=3,5, 
and 7) are due mostly to the variation of the band sepa- 
ration. As shown in the figure the positions of the nodes 
correlate well with the experimental values (open circles) 
taken from the upper panel of FigEJ 

The quantal heating produces an additional node of the 
dc bias induced oscillations. It changes the systematic 
placement of the node positions described by Eq. fp~3|) . In 
the case of a strong quantal heating (as in Fig|9|) the addi- 
tional node occurs at the very beginning of the resistance 
evolution. Expected from Eq. (fT3| node counting can be 
largely restored by a reduction of the node index by two, 
which is the difference between consecutive indexes k in 
Eq. (|13p . The corresponding transformation is shown in 
FigOHl the dashed line is the shift by two units to the 
right of the solid line representing index k in Eq. ([H 
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CONCLUSION 

Quantum oscillations of nonlinear resistance, which oc- 
cur in response to electric current and magnetic field 
applied perpendicular to GaAs quantum wells with two 
populated subbands, are investigated. At small magnetic 
fields the current-induced oscillations are found to be re- 
lated to Landau-Zener transitions between Landau levels 
inside the lowest subband. The period of these oscilla- 
tions is proportional to the magnetic field. At high mag- 
netic fields a different kind of quantum oscillations are 
observed. With respect to the dc bias these resistance 
oscillations are quasi-periodic with a period that is inde- 
pendent of the magnetic field. At a fixed electric current 
the oscillations are periodic in inverse magnetic field. The 
period is independent of the dc bias. The proposed model 
considers these oscillations as a result of joint effect of the 
Joule heating in the systems with discrete spectrum and 
the spatial variations of the energy separation between 
two subbands, which is induced by the electric current. 
Obtained results indicate the feasibility of considerable 
modification of the electron spectrum by applied electric 
current in two dimensional electron systems. 
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